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2 ■ ABSTRACT 

\ In this paper we use the predicted spatial development of helical structures along an 

expanding jet to model observed structures and motions in the 3C 120 jet. New results 
of VLBI imaging of the parsec-scale radio jet in 3C 120 at 5 GHz are examined along 
with older long term monitoring results at 5 GHz and older results obtained at 22 and 
43 GHz. The high frequency observations provide detailed information on motions and 
^ ■ structure from 0.5 to 10 mas from the core and the lower frequency observations from 

1 to 30 mas from the core. Proper motions of helical components associated with the 
pattern and of other components that move through the pattern provide estimates of 
flow and helical pattern speeds. Theoretical modeling of the motion and appearance of 
^ ■ the helical pattern allows determination of sound speeds as a function of the jet viewing 

angle. The jet sound speed declines although probably not as fast as adiabatically. At 
a 12° viewing angle the most likely scenario involves a decline in jet sound speed from 
c/3 < aj < c/\/3 at ~ 0.5 mas from the core to 0.1c < aj < 0.25c at ~ 25 mas from 
the core accompanied by some acceleration in the jet flow from Lorentz factor 7 < 5 to 
7 > 7. The sound speed in the cocoon medium around the jet is less well determined 
but is less than the sound speed in the jet probably by a factor of 1.5 - 5. A largest 
possible viewing angle of 15° implies a jet sound speed at the upper limit of these 
estimates and somewhat higher flow Lorentz factors. However, jet morphology argues 
against viewing angles larger than 12°. At smaller viewing angles the jet sound speed 
is lower and at a 6° viewing angle the jet sound speed is about a factor 2 less but the 
flow Lorentz factor is comparable. The decline in radio intensity is on the order of what 
would be associated with isothermal jet expansion. Knot interknot intensity variations 
are greater than would be expected from adiabatic compressions associated with the 
helical twist and we infer the presence of a shock along the leading edge of the helical 
twist in addition to shock or density structures flowing through the helical pattern. Our 
results imply that the macroscopic heating of the expanding jet fluid is less than the 
microscopic energization of the synchrotron radiating relativistic electrons. 
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Subject headings: galaxies: individual (3C 120) — galaxies: jets — galaxies: active - 
radio continuum: galaxies — hydrodynamics — relativity 

1. Introduction 

The radio source 3C 120 is dominated by a variable core and a prominent one-sided jet ex- 
tending from the core on subparsec scales to about a hundred kiloparsecs (Walker, Benson, & 
Unwin 1987). The galaxy has a redshift z = 0.033 (Baldwin et al. 1980), and with H Q = 
h 100 km s _1 Mpc -1 and h/0.7 = hjo = 1 there are ~ 0.65 pc mas -1 . Recent high frequency 
VLBI observations at 86, 43 and 22 GHz (Gomez et al. 1998; Gomez, Marscher, k Alberdi 1999; 
Gomez et al. 2000, 2001; Marscher et al. 2002) have produced high resolution images of the in- 
nermost features and motions in the jet. The 43 and 22 GHz observations, show that the inner 
jet structure, within 10 mas from the core, is complex and exhibits a subluminal to superluminal 
transition on very small scales. Global VLBI observations at 5 GHz and 10.7 GHz (Walker et al. 
2001) find superluminal motions consistent with those observed at 22 and 43 GHz. The 3C 120 jet 
shows a variety of knots and side-to-side structures, and it has been suggested that the jet structure 
of 3C 120 inside 10 mas may be described by components interacting with the external medium 
in their motion along an underlying helical twist (Gomez et al. 1998, 2000, 2001). VLBI images 
at 1.7 GHz, see Figure 6 in Walker et al. (2001), also show structure suggestive of a helical twist 
beyond 25 mas but with considerably longer wavelength. 

Helical structure in relativistic jets can arise as a result of ordered variation in the flow direction 
at the central engine, e.g., precession, and/or as a result of random perturbations to the jet flow 
such as the jet cloud interaction discussed in Gomez et al. (2000). Initial random perturbations 
can trigger or propagate as pinch, helical, or higher order normal modes of jet distortion. Any 
normal mode structure will depend on the initial excited wave frequency or frequencies and initial 
amplitudes, and on the subsequent propagation and growth or damping of these wave frequencies 
along the jet, e.g., axisymmetric perturbations as in Hardee et al. (1998) and Agudo et al. (2001) or 
some combination including asymmetric perturbations as in Aloy et al. (2003) and Hardee & Hughes 
(2003). In these simulations linear and/or non-linear perturbations were applied. In particular, an 
injection event can coexist and interact with a precession induced helical twist (Aloy et al. 2003). 
In general, non-linearly triggered normal modes are damped until reaching a "linear" regime and 
both non-linear triggers and "linear" normal mode structures can coexist on the jets. Within the 
"linear" regime normal mode structures can be satisfactorily modeled by the linearized relativistic 
hydrodynamic (RHD) equations (Hardee et al. 1998; Agudo et al. 2001; Hardee & Hughes 2003). 
The behavior of normal modes and, in particular, helical twist propagation and growth or damping 
can sensitively depend on jet speed via the Lorentz factor, on the sound speeds in the jet and 
surrounding material, and on the rate of jet expansion (Hardee 2003). 

That the normal wave modes predicted theoretically operate on AGN jets has been suggested 
by Lobanov &: Zensus (2001) who successfully fitted the twisted emission threads observed in the 
3C 273 jet by a combination of helical and elliptical surface and internal normal modes. Other 
similar fitting has been performed by Lobanov, Hardee, & Eilek (2003) in the context of M 87. 
Helical twists have been invoked to explain other observed jet structure and, for example, it has 
been shown that superluminal motions and accelerations along a curved trajectory can be produced 
by helical jet models, e.g., the radio source 3C345 (Hardee 1987; Steffen et al. 1995). While some 
fitting of helically twisted structure to observed jet morphology has been performed, to date no 
detailed self-consistent models have been constructed. 



1 The National Radio Astronomy Observatory is a facility of the National Science Foundation, operated under 
cooperative agreement by Associated Universities, Inc. 
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Only in the past few years have time dependent relativistic hydro dynamical codes become 
readily available (Duncan & Hughes 1994; Marti, Miiller, & Ibahez 1994; Koide, Nishikawa, k, 
Mutel 1996; Falle & Komissarov 1996; Aloy et al. 1999a). More recently this has allowed fully 
3D relativistic jet simulations to be performed with resolution sufficient to compare structure with 
theoretical predictions and/or with features observed in AGN jets (Aloy et al. 1999b, 2000; Hardee 
et al. 2001; Hughes et al. 2002; Hardee & Hughes 2003; Aloy et al. 2003). Ideally we would like to 
model observed jet structures numerically and include non-linear effects. In general, computational 
constraints make it very difficult to simulate and model the super luminal jets that must lie near to 
the line of sight. Satisfactory numerical modeling would require very high grid resolution across the 
jet to reduce the effects of numerical viscosity (~ 100 computational zones) Perucho et al. (2004) 
and additional zones outside the jet, lengthy computational grids in order to account for projection, 
and very large storage capabilities in order to account for light travel time effects, e.g., Gomez et 
al. (1997); Aloy et al. (2003), that make it impractical to conduct numerical modeling involving 
many trials. This difficulty can be overcome using theoretical models based on the linearized fluid 
equations. 

In this paper we consider high frequency VLBI images and proper motions from Gomez et 
al. (1998); Gomez, Marscher, & Alberdi (1999); Gomez et al. (2001) that cover the inner 10 mas 
of the jet, new 5 GHz images that cover the inner jet out to 30 mas, and proper motions from 
Walker et al. (2001) that cover this portion of the jet. With these observations, we begin to 
discern what trajectory components follow, whether motions are constant for a given feature and 
between features, and whether there are similarities between component speeds, trajectories, and 
structure at parsec and tens of parsec scales. We use the observed structure and proper motions 
to constrain the macroscopic properties of the jet and surrounding medium. In §2 the observed 
relevant observations are summarized. In §3 we show how to model relativistically moving helical 
twists along with the accompanying line of sight appearance and model the innermost 10 mas 
portion of the jet as a single frequency helical twist on an isothermally expanding jet. In §4 we 
extend the modeling to larger distances with dual harmonic twist frequencies and also consider 
adiabatic expansion. We conclude in §5 with a summary and discussion of the implications of our 
intensity and dynamical modeling. 

2. Observed Structure and Motion inside 30 mas 

The inner 10 mas structure of 3C 120 has been interpreted as determined by the evolution of 
superluminal components and their complex interactions with the external medium and/or under- 
lying jet (Gomez et al. 1998; Gomez, Marscher, & Alberdi 1999; Gomez et al. 2000, 2001; Marscher 
et al. 2002) The association of some of these components with injection of material into the jet, 
leading to the formation of shocks, has been established by simultaneous radio and X-ray observa- 
tions. Dips in the X-ray light curve are interpreted as caused by the disappearance of a section of 
the inner accretion disk past the event horizon of the black hole, while the remainder of the disk 
material is injected into the jet, leading to a flare in the radio light curves and the appearance of 
a new superluminal component in the jet (Marscher et al. 2002). Other components seem to be 
generated by the passage down the jet of such leading shocks, as has been considered for compo- 
nents m, r, and s (Gomez et al. 2001). These are observed to appear in the wake of the strong 
superluminal component o, and present a similar behavior to the trailing components obtained in 
numerical simulations by Agudo et al. (2001). 

Continued monthly monitoring of the inner structure at 22 and 43 GHz with the VLBA has 
revealed rapid changes in the total and linearly polarized intensity, accompanied by a rotation of 
the magnetic polarization vector, interpreted as resulting from the interaction of the jet components 
with the external medium (Gomez et al. 2000, 2001). In particular, these high frequency observa- 
tions also show evidence for the existence of an underlying helical jet structure, which may also be 
affecting the jet's structural evolution. Component o (Figure 1 at 22 GHz) can be resolved into 
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two subcomponents ol and o2 (Figure 1 at 43 GHz) that show changes in the position angle and 
magnetic polarization suggesting interaction with an underlying helical twist (Gomez et al. 2001). 
In general, there are "northern" and "southern" jet components, see Figures 2 in Gomez et al. 




2 4 6 8 

Relative Right Ascension (mas) 



Fig. 1. — Multiepoch VLBI images of 3C 120 at 43 and 22 GHz with the various components labeled. At 43 GHz 
the contour levels start at 4.4 mjy beam - with increase by factors of 2 from there. The intensity peaks at 0.37 
Jy beam -1 (Nov 1997) and 0.33 Jy beam 1 (Feb 1999). The convolving beam is 0.35 x 0.16 mas extended in position 
angle near 9 degrees. At 22 GHz the contour levels start at 4.2 mjy beam -1 except for Dec 1996 which starts at 
2.1 mjy beam -1 and increase by factors of 2 from there. The intensity peaks at 0.67 Jy beam -1 (Dec 1996), 0.38 
Jy beam -1 (Nov 1997), and 0.43 Jy beam -1 (Feb 1999). The convolving beam is 0.6 x 0.3 mas extended in position 
angle near 9 degrees. 

(2000, 2001), that seem to be ejected with different position angles. The southern components are 
usually brighter than the northern components and this could be the result of the overall curvature 
of the jet to the north, e.g., Walker et al. (2001), leading to higher pressures along the southern 
edge combined with projection effects associated with a helical twist of the jet. 

It is tempting to identify prominent southern components as indicating an overall helical twist. 
The spacing between the prominent southern components h, d, and a (22 GHz) is similar and other 
components can appear to outline a projected helical twist between these prominent southern 
components. For example, note the curved path in intensity contours between components d and a 
(December 1996) and between components h and d (November 1997). That prominent components 
always evolve towards the southern side of the jet suggests that projection plays a significant role 
in their location. These "helical" components with similar spacing and motion would indicate 
the wavelength and motion of the helical twist. Other non-helical components would outline and 
indicate flow through the helically twisted jet. 

Components h, d, and a have separations of 2.5 — 3.5 mas at core distances from 2 — 8 mas 
with a suggestion of increasing separation as core distance increases. Observed proper motions are 
(h) 1.75 and (d) 1.71 mas yr _1 (Gomez, Marscher, &: Alberdi 1999; Gomez et al. 2001), and we use 
an average proper motion of ~ 1.73 mas yr _1 for these components as representative of motion of 
the helical twist. The relatively constant spacing and motion along with a "saturated" component 
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intensity structure beyond 2 mas suggests that the motion and spacing of these components with 
(5^ < 4, where 1 mas yr 1 = 2.2h 7 Q c, should be associated with a saturated helical twist in the 
high frequency regime (Hardee 2003). Here high frequency is relative to the fastest growing or 
"resonant" frequency of the Kelvin-Helmholtz unstable helical twist. 

Non-helical components can move faster than the helical components. For example, in Gomez 
et al. (2001) fast moving components between o and h (e.g., / and 12 in Figure 1) show motions 
from 2.17 to 2.48 mas yr -1 . A wide variety of motions, from ~ 1.5 — 3.5 mas yr -1 are found by 
Walker et al. (2001) at 5 and 10.7 GHz out to 14 mas from the core and are shown along with 
estimated errors in Figure 2. An average motion associated with fast moving components might 
be ~ 2.7 mas yr" 1 . Here this estimate lies between the 2.17 mas yr -1 from Gomez et al. (2001); 
Walker et al. (2001) and the 3.47 mas yr" 1 from Walker et al. (2001) and represents a simple proper 
motion average of components B to J in Figure 2. The resulting apparent superluminal motion, 
(3j b ~ 6, might be assumed to correspond to the flow speed and constrains the viewing angle to less 
than 19°. 

— i 1 1 1 1 1 1 1 1 1 1 — 

A 1.44 ± 0.46 mas/yr G 2.39 ±1.60 

B 3.14 ± 0.12 H 3.48 ± 0.59 4 




1980 1982 1984 1986 1988 



Time of Observation (Yr) 

Fig. 2. — Motions of components seen at 5 and 10.7 GHz in 3C 120 between 1978 and 1988. The components 
are labeled and a straight line with slopes corresponding to the fitted proper motion is drawn through the points. 
The proper motion for each feature along with formal errors is written on the figure. Component labels here do not 
correspond to components in Figs. 1 & 3. For additional details see the Figure 4 caption and Table 3 in Walker et 
al. (2001). 

Interestingly, component o has a proper motion, 1.83 mas yr" 1 , similar to prominent southern 
components, and spacing between o and h is similar to the spacing between h, d, and a. Component 
o appeared in 1998 after a strong outburst in December 1997 (Gomez et al. 2001; Marscher et al. 
2002). Thus, component o should not necessarily be related to the helical components h, d, and 
a assumed formed by projection of the helical pattern. Also component o is not clearly on the 
southern edge of the jet. Thus, we do not now identify o as a helical component. It will be 
interesting to see if in future it behaves like h, d, and a. Such a scenario might occur if a non- 
uniform outflow is precessed and components trigger /merge with the overall helical twist. In this 
case component shocks can lead to a shock along the leading edge of the helical twist. 

Components s, r, and m, interpreted as structures triggered by passage down the jet of a shock 
associated with component o (Gomez et al. 2001), have motions « 0.27, 0.40, and 0.49 mas yr" 1 
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respectively, and the separation of these components increases from s to o. Acceleration and 
increasing spacing occurs naturally for pinch body mode components triggered near to the resonant 
frequency as the pinches transition from the resonant to the high frequency regime by virtue of jet 
expansion, e.g., Agudo et al. (2001). Alternatively these components could be associated with the 
helical twist. Here again acceleration and increasing spacing results from a change in helical wave 
speed as the twist transitions from the low or resonant frequency regime to the high frequency regime 
by virtue of jet expansion (Hardee 2003). In either case the proper motion of the component nearest 
to the core with ~ 0.6 could represent motion of pinching or helical twisting in the resonant or 
low frequency regime. We note here, and address this issue in §3.1, that the proper motion of pinch 
body mode components or "helical " components associated with the helical surface (twist) mode 
are similar in the resonant and high frequency regimes. Thus, the identification of these innermost 
components either with pinches or helical twist can give a similar final modeling result. 



T 
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Fig. 3. — VLBI images of 3C 120 at 5 GHz from 1999.11 and 1999.72 September 21. The contour levels start with 
-1, -0.5, 0.5, and 1 mjy beam -1 with increase by factors of 2 from there. The peak flux densities are 1.59 and 1.00 
Jy beam -1 in the first and second images respectively. The convolving beams are 3.5 x 1.0 mas extended in position 
angle near 9 degrees. The marked jet features in the 1999.11 image correspond to features seen in the 1999 February 
image of Figure 1, which is based on 22 GHz data from the preceeding day. The core position marked here is at the 
position of the bright, easternmost feature of the 22 GHz image. Clearly the core is significantly absorbed at 5 GHz. 

Inspection of images in Gomez et al. (2000, 2001) shows that the jet is bent to the north 
inside 10 mas with a position angle change on the order of A°p A ~ 12°, although most of this 
change appears to occur beyond 2 mas from the core. At an assumed viewing angle of 9 ~ 10°, 
this apparent change coresponds to a intrinsic change of Apa ~ sin 9 A°p A ~ 2°. The apparent 
and intrisic bending angles corresponding to this position angle change are somewhat larger, i.e., 
the change in angle, A, defined by a vector tangent to the jet is greater than the change in 
position angle, Apa, defined by the vector from the core to a point on the curved jet. An intrinsic 
bending angle of say A = 3°, even if all bending were to occur in the innermost 2 mas could 
not account for the observed acceleration. Provided A << 9 « 1 radian, the change in motion 
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5(3 /(3 ~ A/0 << 1 and jet bending to the north cannot account for the observed acceleration 
from (3 ob ~ 0.6, component s, to (3 ob ~ 1.1, component m, and falls far short of an acceleration to 
(3 oh ~ 4, component o. 

Jet structure out to ~ 30 mas is revealed by the total intensity maps in Figure 3. The 
observational details for these previously unpublished images are given in Appendix A. Here we 
identify the first four components with o, h, d, and a. In Figure 3 intensity contours between h, d, 
and a indicate a helical morphology like that suggested by Gomez et al. (1998) and again with a 
suggestion of increasing spacing (typical spacing ~ 3.5 mas) as core distance increases. Between 
10 and 15 mas there is a discontinuity in the regular pattern of components. If we interpret the 
outermost structure in terms of a helical pattern then the wavelength has shown an abrupt increase 
to ~ 8 mas. This increase in helical component spacing is also accompanied by an apparent increase 
in the proper motions of components where component motions beyond about 15 mas from the 
core determined from 1.7 GHz multiepoch images are between 2.5 and 3.1 mas yr _1 (Walker et al. 
2001). This would imply an increase in observed helical wave speed from (3 ob < 4 to (3 ob ~ 6. 

A 1.7 GHz image shown in Figure 4 shows structure between 30 - 80 mas that is qualitatively 
similar in appearance to the helical morphology inside 30 mas. Note the prominent components 
are again on the southern side of the jet and note a curved path in intensity contours between 50 
- 80 mas from the core. Here we see an indication of a component spacing jump from ~ 15 mas 
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Fig. 4. — A 1997.70 image of the inner 360 mas of the jet at 1.7 GHz. The image has been convolved to a beam of 
12.5 mas x 4.0 mas elongated in position angle -10°. The lower contour levels are -1.5, -0.75, 0.75, and 1.5 mjy beam -1 
with the higher levels increasing by factors of 2 from there. 

between 35 - 50 mas from the core to ~ 30 mas between 50 - 80 mas from the core. Beyond the 
slowly moving component, (3 ob < 1, at ~ 80 mas, the jet character changes significantly. This image 
and unpublished 5 GHz images suggest spacing jumps by factors of two as core distance increases. 
At least approximately helical component spacing makes a transition from an apparent < 3.5 mas 
spacing within 10 mas from the core to ~ 30 mas by 80 mas from the core. The implied increase in 
helical component spacing, approximately proportional to the jet radius, would require a systematic 
linear increase in flow speed as a function of core distance if it was a flow speed effect. On the other 
hand, this increasing component spacing may be representative of excitation of multiple harmonic 
helical frequencies and with lower frequencies supplanting higher frequencies at larger core distance. 

It is from the analysis of the wavelength and motion of helical structure along with the motion 
associated with the jet flow that we expect to constrain, through our modeling, the macroscopic 
properties of the jet and surrounding medium. 
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3. Modeling Helical Structures 



Comparison between theoretical predictions and numerical simulations indicates that helical 
structures must be operating near to the linear regime if jets are to remain highly collimated. In 
the non-linear regime helical jet surface distortion at resonant or longer wavelengths leads to mass 
entrainment that can disrupt highly collimated flow and can lead to the formation of an extended 
velocity shear layer. Note however, that helical surface distortion at wavelengths much shorter than 
resonant may be able to reach a non- linear saturation without disrupting highly collimated flow, e.g., 
Xu, Hardee, & Stone (2000). Even when formally weakly non- linear, normal mode structures can 
be satisfactoriy modeled using the linearized fluid equations (Hardee et al. 1998; Hardee & Hughes 
2003). In the linear or weakly non-linear regime moving helical structures and accompanying 
flow fields found from the linearized RHD equations behave in predictable ways and the resulting 
behavior can lead to estimates of sound speeds from observed proper motions. In what follows we 
show in detail how to use the observed proper motions of features in the inner 30 mas of the 3C 120 
jet to obtain the first order consistent flow field associated with the helical structure, produce 
pseudo-synchrotron intensity images appropriate to adiabatic compressions associated with helical 
structure, and most importantly constrain the parameter regime in which the jet must operate. 



The linearized RHD equations show that a helical twist is Kelvin-Helmholtz unstable and 
propagates at a wave speed dependent on jet speed, sound speeds, and on the frequency of the 
wave relative to a "resonant" or maximally unstable frequency, uj* . For sufficiently supersonic flow, 
i.e., Mj = Vj/aj >> 1 and M x = Vj/a x >> 1, the resonant frequency is uj* ~ 1.5a x /Rj where Rj 
is the jet radius and a^ x is the sound speed in jet or external medium. In general, solution of the 
wave dispersion relation is necessary to compute the wave speed and wavelength when the wave 
frequency is within an order of magnitude of uj* . Along an expanding jet uj* oc RJ 1 decreases and 
solution of the dispersion relation is needed to follow the development of a helical wave of constant 
frequency as it is advected along the jet from the low to the high frequency regime. 

At frequencies more than an order of magnitude below uj* (uj « uj*) the wave speed is given 



a = [rP /(po + [r/(r - l)]P /c 2 )] 1/2 , where 4/3 < V < 5/3 is the adiabatic index. The density, 
po, and pressure, Po, are measured in the proper fluid frames, and since pressure balance has been 
assumed i] is an enthalpy ratio if the adiabatic indices are the same inside and outside the jet. 

The apparent helical wavelength in the observer frame is related to the intrinsic wavelength 
by \ ob = (f3^/f3 w )\ where (3^ = [sin 0/(1 - (3 W cos 0)](3 W is the apparent wave speed, and the 
observed wavelength can be greater or less than the intrinsic wavelength depending on wave speed 
and viewing angle 9. In general, change in the observed wavelength is amplified relative to change 
in the intrinsic wavelength with 



3.1. Moving Helical Patterns 



v w = {w/k)\ Rea i ~ — — y-Vj , 

1 + 7 Tj 

and at high frequencies when uj » uj* the wave speed is given by 




(1) 




\°\z ob ) = X 



(zsm.9) = 



sin 9 



\{z) . 



(3) 



1 — Pw(z) cos 9 
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In equation (3) z is the intrinsic distance to the core where z b = zsinO is the observed distance. 
If the wave speed is relatively low at low frequencies, i.e., 7 2 r/ < 1, then the speed increases until 
it reaches the value given by equation (2) at high frequencies. In this case a helical wave advected 
along an expanding jet can appear to increase in wavelength approximately proportional to Rj as 
the wave speed increases. If 7 2 r/ >> 1 there is a minimum in the wave speed near to the resonant 
frequency and there can still be significant change in the observed wavelength. 

Equation (1) can be rewritten to give the sound speed ratio in terms of the observed flow 
speed, wave speed {uj « uj*) and viewing angle as 



[1- (pf sin6>-cosfl) 2 ]/3f 
' ' ' ~ (sin 9 + Pf cos 9) (f3f - p ob ) sin 9 



V - Ma^ = , . • „ , '"" n • (4) 



Equation (2) can be rewritten to give the jet sound speed in terms of the observed flow speed, wave 
speed (uj » uj*) and viewing angle as 



Pf-PZ b 

(Pf + P$) cos 9 - (PfP ob - 1) sin 6 



P» ~ a i/° ' Rob , Rob\ nna a _ (Rob Rob i ^ ----- n ' ( 5 ) 



Equation (3) can be rewritten to give the intrinsic wavelength in terms of the observed wave speed, 
wavelength and viewing angle as 

In equations (4 - 6) fi oh = v ob /c is the apparent super or subluminal flow or wave speed, and we 
have used p = [p ob /(sm9 + p ob cos9)]. Provided pfp ob > 1, the superluminal case considered here, 
a lower limit to the jet sound speed 



Pf-Pt 
Pf + Pt 



can be found from equation (5) for sm.9 — > 0. As the viewing angle increases the sound speed 
increases and in the limit P s — > l/\/3 a maximum viewing angle 

^ _ ^3-1 

Vmax - ^ {&) 

is obtained from equation (5) where we have used cos# ~ 1, sin# ~ 9, and assumed pfp ob » 1. 

It can be shown that the supersonic limit, PJP f < l/y/3, guarantees that P° b I Pf > (1 - P s /Pf) 

and 9 max > 0. As P ob -> Pf the less restrictive 9 max ~ 2/ 1 P° b is recovered from equation (8). 

Given that the observed prominent "helical" component proper motions in the inner 3C 120 
jet are well defined, the largest uncertainty in the modeling process lies in the assumed flow speed 
which could be as small as P ob | m i ra ~ 5 from the observed proper motions of ~ 2.3 mas yr _1 found 

between components o and h at 43 GHz (Gomez et al. 2001) or as large as P° b \ max^ 7.5 from the 

largest observed proper motion of ~ 3.4 mas yr -1 at 5 GHz (Walker et al. 2001). The section of 
the jet that we consider also clearly contains regions of different component spacing and motion. 
Thus, we are faced with the possibility of flow acceleration and changing sound speeds along the jet. 
Initially we shall consider two cases. In both cases let us assume P ob (u) << uj*) ~ 0.6 but otherwise 
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Pw(uj » co*) ~ 4 (appropriate to inside 10 mas) or ^{uj » 



obi 



5.5 (appropriate to outside 



10 mas). Here we are identifying the subluminal motions of the innermost observed components in 
the 3C 120 jet with a minimum helical wave speed. Equations (4) and (5) then provide preliminary 
estimates of the sound speeds and Figure 5 shows how j3 s = aj/c and a x /c = rj l / 2 [5 s vary as a 
function of viewing angle, 9, and observed flow speed, (3 ob . Examination of Figure 5 indicates that 
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Sound speeds from eqs.(4) & (5) as a function of viewing angle and apparent jet flow speed f]f b = 5 
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Fig. 5 

Jet sound speed for two high frequency wave speeds, (w >> w*) ~ 4 (left top panel) & 5.5 (right top panel). The 
value of the external sound speed that accompanies the high frequency wave speed is shown in the panels immediately 
below. Here in both cases we assume that the observed low frequency wave speed is ffi (o> << w*) ~ 0.6. 

the external sound speed would lie within the range 0.05 < a x /c < 0.15 and the internal sound 
speed aj/c > 0.12 for ^(uj » u*) ~ 4. Values are 0.02 < a x /c < 0.08 and aj/c > 0.04 or factors 



2 - 3 times lower for @%(u » oj*) ~ 5.5. Equations (7) and (8) with ft. 



ob 



4 and pf 



6 suggest 



that the viewing angle is 9 < 15° and we will assume that 9 = 15° represents the largest possible 
viewing angle. The fact that implied sound speeds are lower for the higher assumed wave speed 
beyond 10 mas suggests a jet cooled by expansion. 

When the observed spatial change in wave motion and wavelength imply a helical wave with 
frequency within an order of magnitude of resonance, the full wave dispersion relation must be 
used to fit the observed wave motion and wavelength, and to estimate sound speeds. In our initial 
modeling we will follow the assumption that the jet can be modeled as an isothermal constant speed 
expansion with isothermal external medium. For viewing angles 9 = 15°, 12°, 9°, 6°, and 3° we have 
considered two different cases that bracket the likely flow and sound speed range. In both cases 
we consider the portion of the jet where the jet radius ranges from 0.046 mas < Rj < 1.56 mas 
(0.03 pc < Rj < 1.01 pc), the observed distance ranges from 0.44 mas < < 14.88 mas, and the 
observed half-opening angle ip Q b = 6° w 0.105 radian. 
For Case 1 we assume a flow speed j3j b « 6 with: 

• The observed wave speed « 0.7 at z b = 0.44 mas and /3° 6 

• At z b = 14.88 mas, the observed helical wavelength is X ob 

Here the observed wave speed ranges from typical of the proper motions of the innermost compo- 
nents s ((3 ob ~ 0.6 @ z a b ~ 0.33 mas) and r (f3 ob ~ 0.9 @ z a b ~ 0.55 mas), but does not achieve the 
proper motion of the "helical" components (/3° 6 ~ 4) for z Q b ~ 2 — 10 mas or the proper motion of 
the "helical" components for z Q b > 15 mas (/3° 6 ~ 6). Additionally, the predicted wavelength for 



4 at z ob = 14i 
3.6 mas. 



mas. 



2.3Rj ~ 
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z b ~ 2 — 10 mas is somewhat less than the observed wavelength. At least approximately this case 
is representative of the jet close to the core. 

For Case 2 we assume a flow speed (3°^ « 6.5 with: 

• The observed wave speed /3° 6 « 0.8 at z G f, = 0.44 mas and « 5 at z Q b = 14.88 mas. 

• At z b = 14.88 mas, the observed helical wavelength is X ob ~ 3.1i?, ~ 4.3 mas. 

Here the observed wave speed must begin somewhat above that of the typical proper motions of the 
innermost components {(3^ ~ 0.7) in order to achieve the typical proper motion of the "helical" 
components (/3° b ~ 4) for z b ~ 2 — 10 mas, and to nearly achieve the faster proper motion of 
"helical" components ~ 6) for z^ > 15 mas. Here the predicted wavelength is somewhat 
larger than the observed wavelength for z f, ~ 2 — 10 mas. 

The above set of constraints along with the preliminary estimates of allowed sound speeds, 
e.g., Figure 5, gets us into the correct region of parameter space. In general, we first obtain an 
estimate for the jet sound speed by fitting the wavelength and wave speed at z Q b = 14.88 mas. 
Here, where the wavelength is in the short wavelength and high frequency regime, the wave speed 
is nearly independent of the external sound speed and primarily depends on the jet speed and the 
jet sound speed. Typically this requires 10 or more trials at a single viewing angle. Note that we 
also obtain ujRj/u at this location and with Rj and u known we obtain the angular frequency to. 
Subsequently we find the external sound speed that would yield the wave speed at z a b = 0.44 mas 
and u)R$/u given our estimate of the jet sound speed. Here the angular frequency, uj, of the wave 
is conserved but ujRq/u is reduced relative to ojRj/u by the smaller jet radius. A few additional 
trials may be necessary as the value for the external sound speed modifies slightly the required jet 
sound speed and value of uiRj/u. The resulting estimates for the sound speeds will be accurate on 
the 10% level, i.e., 10% differences provide a similar satisfactory fit. Computing time is minimal, 
basically solving the dispersion relation for different values of the sound speed, but the non-linear 
nature of the dispersion relation near to resonance makes this a procedure that would be difficult 
to automate. 

Figure 6 shows a sample of the dispersion relation solutions that best satisfy the set of condi- 
tions above along with the intrinsic and observed wavelengths and wave speeds and Table 1 contains 
the associated parameters. 

Table 1: Jet Parameters for Case 1 & Case 2 



e 


ip (rad) 


1: vj/c 


7 


CLj/c 


a x /c 


2: vj/c 


7 


dj/c 


a x /c 


15° 


0.0271 


0.9910 


7.47 


0.460 


0.095 


0.9943 


9.37 


0.400 


0.040 


12° 


0.0218 


0.9874 


6.31 


0.330 


0.110 


0.9900 


7.09 


0.250 


0.045 


9° 


0.0164 


0.9864 


6.09 


0.260 


0.130 


0.9885 


6.61 


0.180 


0.050 


6° 


0.0109 


0.9882 


6.53 


0.210 


0.150 


0.9895 


6.92 


0.130 


0.055 


3° 


0.0055 


0.9927 


8.29 


0.150 


0.170 


0.9935 


8.78 


0.090 


0.065 



For Case 1 the observed wave speed, /3° ~ 0.7, corresponds to the wave speed at frequency 
ujRo/vj ~ 0.32. The implied angular frequency is uj ~ 1 x 10~ 7 radian/sec for a periodicity of about 
2 years. The jet sound speed is lower than would be suggested by eq. (5) in order to achieve srf 4 
at z b = 14.88 mas where toRj/vj ~ 11 as the jet expands by a factor ~ 34 over the modeled region. 
The frequency is above "resonance" and thus the value for the external sound speed, found from 
exact solution to the dispersion relation, varies from about 10% less to 40% more than is found 
using eqs. (4) & (5) as viewing angle decreases from 15° to 6°. We note that the indicated external 
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sound speed at the 3° viewing angle would need to be a x /c ~ 0.155 to give /3° 6 | m m~ 0.6. The fact 
that we cannot obtain an external sound speed estimate consistent with the required subluminal 
motions at 9 < 3° and the observed subluminal acceleration would suggest that 6 > 3°. However, 
both observational errors and uncertainties in the interpretation of features makes a firm lower 
limit to the viewing angle unlikely. 
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Fig. 6. — Solutions fci?j as a function of uiRj/vj for helical waves using parameters appropriate to viewing angles 
9 — 15°, 12°, 9° & 6° for Case 1 (top panels) and Case 2 (bottom panels). Dotted (Dash-dot) lines indicate the real, 
k r , (imaginary, ki) part of the wavenumber. The intrinsic (dotted) and apparent (dashed) wave speeds and intrinsic 
(short dash-dot) and apparent (long dash-dot) wavelengths are shown immediately below the appropriate dispersion 
relation solutions. The long, medium and short vertical lines mark ui* , uii(z b ~ 0.45 mas), and uif(z b ~ 14.9 mas), 
respectively. 

For Case 2 the observed wave speed, /3° 6 ~ 0.8, corresponds to the wave speed at frequency 
uRo/vj ~ 0.30. The implied angular frequency remains w~lx 10~ 7 radian/sec for a periodicity 
of about 2 years. The jet sound speed is slightly lower than would be suggested by eq. (5) in order 
to achieve « 5 at z b = 14.88 mas where uRj/vj ~ 10. The frequency is above "resonance" 
and and thus the value for the external sound speed, found from exact solution to the dispersion 
relation, varies from about 30% less to 15% more than is found using eqs. (4) & (5) as viewing 



angle decreases from 15° to 6°. For this parameter set (5, 



ob 



t < 0.6 can be obtained at all viewing 
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angles and no lower limit to the viewing angle is implied by the observed subluminal acceleration. 

In the above we have assumed that the subluminal motions of the innermost 3C 120 jet com- 
ponents would be representative of the helical wave speed. This would not be true if the helical 
wave was in the low frequency limit at this location because low frequency propagation of helical 
surface and pinch body modes can be very different. However, we find that the helical wave at this 
location is somewhat above resonance. A normal mode analysis shows that wave propagation at or 
above the resonant frequency for pinch body modes and for helical twist modes is approximately 
the same (Hardee 2000) . Since our estimate is based on a helical wave with initial frequency above 
resonance, there will be only modest difference in the estimated external sound speed if the sublu- 
minal component motions are modeled as pinch body mode structures. The variation in external 
sound speed from this initial position has little influence on our results. This is because the modeled 
helical structure primarily lies in the high frequency regime where wave speed is nearly independent 
of the external sound speed or its variation. 

3.2. Modeling Helically Twisted Fluid Flow 

As a helical twist is Kelvin-Helmholtz unstable, the displacement amplitude of the jet surface 
grows according to 

A = A exp / l(z)~ x dz , (9) 

where Aq is the displacement amplitude at zq and l{z) = k^ 1 (z) is the spatial growth length. 
Growth at the predicted linear growth rates has been verified in the "linear" regime by non- 
relativistic simulations (Stone, Xu, & Hardee 2000; Xu, Hardee, & Stone 2000) and by relativistic 
simulations (Perucho et al. 2004). The magnitude of accompanying velocity and pressure fluctua- 
tions grows according to A(z)/Rj(z), and for constant jet expansion 



A(z) A ex P 



(10) 



Rj(z) R [l + (z-zoW/Ro] ' 
where Rj = Rq + (z — zo)ip. For ui < uj* the growth length is given by 

£(u) = 'yr} 1/2 (uR j /v j )~ 1 R j = jMjiuRj/a^Rj . (11a) 
At uj ~ uj* the minimum growth length 

t > ^M.R, (lib) 

is somewhat larger than a straight extrapolation of the low frequency result to the resonant fre- 
quency. At frequencies 1 < uj/uj* < 10 the growth length rapidly increases with 

t{u) ~t(uj/uj*) m , (11c) 

where m < 1. At higher frequencies the growth rate can plateau depending on jet and sound 
speeds. The dependence of the growth length on wave frequency relative to uj* means that a 
maximum in A(z)/Rj(z) occurs on an expanding jet when uj > uj* . Note that the jet half opening 
angle, tp, must be less than the relativistic Mach angle, ~ (7Mj) _1 , for self-consistency. Equation 
(10) reflects this self-consistency in that A(z)/Rj(z) < Aq/Rq as ip — > (^Mj)" 1 and no growth in 
pressure or velocity fluctuation is predicted. 

Figure 7 shows how pressure and velocity fluctuations grow for a viewing angle of 12°. The 
ID cuts are made at locations 0.22 < x/Rj < 0.88 and at y = 0. Thus, v x is a radial velocity 
and v y is a toroidal velocity component. Calculation of the pressure and velocity structure has 
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Fig. 7. — The topmost panel of each set shows the growth or decline of A/Rj as a function of distance along the 
jet axis for a viewing angle of 9 = 12°. Three panels below show the normalized pressure P/Po(z), and velocity 
components "/ z v z /c, v x /c and v y /c. The vertical dotted lines mark the projected separation from the core in mas. 
ID cuts are at x/Rj = 0.22 (dash), 0.44 (long dash), 0.66 (short-long dash), 0.88 (solid). 

used expressions given in Hardee (2000). The spread in radial velocity from jet center to jet 
surface indicates jet expansion and there is considerably more fluctuation near jet center than 
at the jet surface. Uniform v y across the jet indicates uniform sideways motion of the jet. The 
axial velocity shows more fluctuation at the jet surface than near the jet center. For all angles 
we have allowed the growth of perturbations from some initial amplitude Aq/Rq at position zq 
up to a saturation amplitude assumed equal to the maximum allowed by the linear analysis, i.e., 
P ~ Po(z) ± Po(z). For all viewing angles and for both parameter sets we find that an initial 
amplitude can be chosen such that saturation is achieved without imposing an arbitrary decline in 
the growth rate at high frequency. Thus, we find that saturation can occur naturally as a result 
of decline in the growth rate at high frequency combined with jet expansion at constant opening 
angle. For Case 1 parameters where ujRo/vj > iv*Ro/vj, A/Rj increases before declining but for 
Case 2 parameters where cuRo/vj >> uj*Ro/vj, A/Rj declines. Nevertheless, the initial velocity 
and pressure perturbation at z Q b < 0.4 mas that can lead to saturation is very small. In general, for 



-15- 



a given pressure fluctuation, velocity fluctuations are reduced for lower values of the jet sound speed 
and higher values of the Lorentz factor, i.e., higher relativistic Mach number 7M7. A saturation 
amplitude of AP ~ Po corresponds to flow with helical pitch at about half the relativistic Mach 
angle. For saturated helical waves with oj » lj* displacement of the jet surface A << OARj. 

3.3. Modeling Helical Intensity Morphology 

In intensity modeling we do not attempt to reproduce details within the jet or the exact 
positioning of "helical" components. Line of sight images are constructed assuming that a pseudo- 
synchrotron emissivity at fixed frequency can be written as 

e„ oc n)- 2a p] a (B sin 9 B ) 1+a D 2+a (12) 

where 9b is the angle of the magnetic field to the line of sight and D = [7(1 — (3 cos 6?)] _1 is the 
Doppler boost factor for fluid flowing with speed (5 at angle 9 to the line of sight. The theoretically 

computed flow fields are used to achieve correct inclusion of Doppler boosting effects on intensity 

2 /3 

variations. The magnetic field strength is assumed to decline and fluctuate as B oc n • appropriate 
to a disordered magnetic field. This form for the pseudo-synchrotron emissivity is one way to 
calculate intensities when the synchrotron emitting particles are not explicitly tracked (Clarke, 
Norman, & Burns 1989), and when radiative losses and acceleration processes are not included. 
Jones, Ryu, & Engel (1999) have shown that pseudo-synchrotron intensities computed this way are 
not too different from synchrotron intensities that include radiative losses and shock acceleration 
processes, and so can provide an acceptable image of jet structure. In our intensity modeling we 
assume that the magnetic pressure is negligible and the jet expands at fixed opening angle. In what 
follows we will model the jet assuming isothermal constant speed expansion, i.e., conserve energy 
flux, and that the sound speed ratio, a x /aj is constant. Here the jet Mach number, Mj = vj/aj, 
remains constant and the external Mach number, M x = Vj/a x , also remains constant with the jet in 
pressure balance with an external cocoon medium. Along the conical constant velocity isothermal 
jet the particle number density and pressure decline oc RJ 2 - For the 3C 120 jet we use a ~ 0.65 
where I{y) oc v~ a . 

Figure 8 shows the resulting pseudo-synchrotron intensity images for four viewing angles. Here 
the intrinsic jet becomes longer as the viewing angle decreases so that the images represent the 
observed projected jet. All light travel time effects and Doppler boosting effects are included. The 
brightest radio emission traces the path of the helically twisted high pressure ridge which effectively 
lies within the conical jet's surface. The helical pitch of the high pressure ridge is decoupled from 
the helical pitch of the flow primarily as a result of rapid wave motion when u » lo*, and the 
intrinsic flow pitch angle is less than would be inferred from the intrinsic helical pitch of the high 
pressure ridge. Line of sight effects mean that the brightest regions will be on the south (lower) 
side of the jet where the high pressure ridge approaches closer to the line of sight and increased 
somewhat by increased Doppler boosting. On the north (upper) side of the jet the pressure ridge 
is stretched by projection and the brightness is decreased somewhat because of reduced Doppler 
boosting. The similarity in appearance for the different viewing angles occurs because the flow field 
and Doppler boost factor exhibit less angular variation at the smaller viewing angles. This reduced 
angular flow variation is a result of the lower jet sound speed and higher Lorentz factor required at 
the smaller viewing angles for the observed flow speed. 

Helical component formation is enhanced at the shorter wavelengths of Case 1 and is a combi- 
nation of projection and enhanced Doppler boosting. However, our intensity variations which are 
the result of projected adiabatic compressions show knot interknot variation less than a factor of 
4. This is less than the variation indicated by the 5 GHz images (Figure 3), and considerably less 
variation than indicated by the 22 GHz and 43 GHz images (Figure 1). This result suggests that 
compressions resulting in the "helical" component formation in the 3C 120 jet are not adiabatic. A 
secondary brightening on the northern side of the jet between the brighter southern components is 
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Fig. 8. — Case 1 (left column) and Case 2 (right column) pseudo-synchrotron intensities for viewing angles 15°, 
12°, 9° and 6° (top to bottom). Contours are in factors of y/2. Gray scale relates to intensity but gray scale and 
contours are adjusted relative to the maximum intensity in each image and cannot be intercompared. 

a result of projection effects and relativistic aberration associated with motion of the ribbon like 
high pressure region. At the smaller viewing angles the flow angle relative to the line of sight is 
sufficiently increased at this location to reduce Doppler boosting and reduce this secondary peak. 
The fact that the 22 GHz intensity structure indicates a secondary brightening on the northern 
side between helical components suggests that the 3C 120 jet lies at larger viewing angles. 

4. Multiple Helical Frequencies and Adiabatic Expansion 

The jump in component spacing seen in the 5 GHz images (Figure 4) cannot be reproduced by 
a simple increase in observed wavelength of a helical twist of single frequency. We have investigated 
a combination of harmonic frequencies, LV2 ~ 2wi, for the isothermal jet expansion model and also 
for an adiabatic jet expansion model in which flow acceleration occurs. In what follows we restrict 
our attention to the single viewing angle of 12°. 

For a single frequency excited by precession of the central engine the behavior of a helical twist 
is predictable over many orders of magnitude of jet expansion. Difficulties arise if multiple frequen- 
cies are excited. The largest uncertainty in dealing with multiple frequencies lies in the interaction 
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between multiple frequencies associated with the same wave mode. Numerical simulations have 
shown that the faster growing higher order normal modes (elliptical, triangular etc.) saturate, do 
not slow the growth of a helical mode wave, and decline in amplitude as the longer wavelength 
helical twist grows (Hardee, Clarke & Rosen 1997). Another simulation shows that pinch mode 
waves show abrupt change to longer wavelength on a conically expanding axisymmetric jet (Agudo 
et al. 2001). Therefore we follow the assumption that a high frequency saturated wave will decline 
in amplitude as a lower frequency wave grows as specified by the linear theory. 

4.1. Isothermal Expansion Model 

We use the isothermal jet parameters of Case 2 as we are interested in structure beyond 10 mas 
and consider uiRq/vj ~ 0.18 and U2Ro/vj ~ 0.36 at ~ 0.44 mas. This choice of frequencies 
provides the closest approximation to the wavelengths evident in 43, 22, and 5 GHz images (Figs. 1 
& 4). The initial wave amplitudes, A2/R0 ~ 0.04 and A\/Rq ~ 0.03, have been chosen so that the 
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Fig. 9. — The topmost panel of each set shows the growth or decline of A/Rj as a function of distance along the 
jet axis for a viewing angle of 9 = 12°. The dotted lines & dashed lines show A^jRj & Ai/Rj, respectively. Three 
panels below show the normalized pressure P/Po(z), and velocity components ^f z v z /c, v x /c and v y /c. The vertical 
dotted lines mark the projected separation from the core in mas. ID cuts are at x/Rj = 0.22 (dash), 0.44(long dash), 
0.66 (short-long dash), 0.88 (solid). 
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lower frequency grows without constraint and as specified by the linear theory. The high and low 
frequency waves are allowed to grow at small core distances as specified by the linear theory until 
the combined amplitudes are saturated. At larger distances the high frequency wave is damped 
as the low frequency wave continues to grow as specified by the linear theory in order to keep 
pressure fluctuations within the "linear" limit. This damping is ad hoc as the linear equations 
used here provide no mechanism for damping. The wave amplitudes along with the accompanying 
pressure and velocity fluctuations are shown in Figure 9. In this self-consistent calculation the 
initial low frequency amplitude achieves pressure fluctuations of about 80% of saturation at large 
core distances. 

A pseudo-synchrotron intensity image for the isothermally expanding jet that corresponds 
to the fluctuations shown in Figure 9 is shown in Figure 10. This image shows modest helical 
components at the longer wavelength, A o6 (u>i) ~ 8 mas when z Q b > 11 mas, but with the less 
than maximal short wavelength amplitude the modest helical components that are evident when 
z b < 8 mas in Figure 8 do not appear. In part the lack of significant helical components is a 
result of both this choice of viewing angle and the lower jet sound speed required by Case 2, i.e., 
higher relativistic Mach number reduces differential Doppler boosting effects. The higher frequency 
wavelength shows steady wavelength increase and X ob (oj2) ~ 2 — 4 mas when z t, ~ 1 — 11 mas. It 
is clear that the basic 3C 120 jet helical structure and motion can be reproduced in this fashion 
and a larger or smaller viewing angle can recover modest helical components in the inner 10 mas. 
However, on the 3C 120 jet the higher frequency must saturate close to the core while still allowing 
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Fig. 10. — Pseudo-synchrotron intensity for an isothermally expanding jet with harmonic frequencies at viewing 
angle 12° corresponding to the fluctuations shown in Figure 9. Contours are in factors of 2. 

the lower frequency to saturate farther from the core as the higher frequency is damped. This 
requires that the higher frequency have a larger initial growth rate than the lower frequency and 
this is not the case for an isothermal expansion model with U2 > oj\ > oo* and ki(i02) < ki(ui) at 
z b > 0.44 mas. Of course, the correct ordering of growth rates will be the situation closer to the 
core. Nevertheless, the high frequency wave must be damped too close to the core and the low 
frequency wave does not reach saturation for our best fit isothermal expansion model. 

4.2. Adiabatic Expansion Model 

The rapid initial growth in wave amplitude for the Case 1 isothermal expansion model assumed 
to represent conditions nearer to the core and the higher wave speeds for the Case 2 isothermal 
expansion model assumed to represent conditions farther from the core suggests that a combination 
of jet sound speed decline and jet acceleration could reproduce the observed wavelength and proper 
motion behavior of helical components along the inner 30 mas of the 3C 120 jet. Here we might 
assume that the conditions associated with Case 1 are representative of the innermost 2 mas of 
the observed jet. To test this hypothesis we have constructed a self-consistent adiabatic expansion 
model for a viewing angle of 12° with parameters varying as indicated in Figure 11. 
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Fig. 11. — The top panel shows the decline in jet density rij(z)/nj(zo) and jet pressure Pj(z)/Pj(zo) oc rij(z) 13 ^ 9 
along with (Ro/Rj) 2 for reference. The bottom panel shows the Lorentz factor, ■y(z), the jet sound speed, [3 S — dj(z)/c, 
the observed flow speed, f3° b (z), and the observed helical component speed in the high frequency limit, f3^ b >>ul *(z). 
The vertical dotted lines mark the projected separation from the core in mas. 

To produce this model we have used the adiabatic wind equation 



1 dj 



along with the continuity equation 



and 



a(z) 



7f3nz 2 = <y (3 n z 2 > 



p{z)c 2 + _r_p (z) 



1/2 



(13) 
(14) 

(15) 



where P (z) oc n(z) r with T = 13/9 appropriate to a mixture of hot electrons and colder baryons 
(Synge 1957), and P /p c 2 = 0.28 to give a /c = 0.460 at z ob = 0.44 mas and with vq/c = 0.9650. 
In the above, constant jet expansion at fixed opening angle is assumed, i.e., Rj oc z. 

Along the adiabatically expanding jet the growth and wavelength of helical frequencies must 
be tracked explicitly as conditions in the jet and in the external medium change. In calculations 
the sound speed ratio a x /aj = 0.75, which gives ~ 0.7 at z \, = 0.44 mas, will be assumed to 
remain fixed as dj declines. There is no particular physical reason for this choice but results are 
only weakly dependent on this assumption. An illustration of the growth rate, observed wave speed 
and observed wavelength for three different initial harmonic frequency pairs, uj<i ~ 2uj\ is shown in 
Figure 12. All cases begin at z Q b = 0.44 mas (zq ~ 1.4 pc) where i?o ~ 0.03 pc and vq/c = 0.9650. 
The position z \, ~ 25 mas (z ~ 80 pc) is indicated by vertical lines in the figure and corresponds to 
ujRj/vj ~ 55ujRq/vo. The three panels span the possible frequency range that will reproduce the 
observed wave speed and wavelength behavior between the initial and final position. Case C begins 
with the high frequency UJ2 = ^* a t -z fc = 0.44 mas and is approximately the highest frequency for 
which ki{uj2) > ki{uj\). In particular, case B with loiRq/vq = 0.214 and c^-Ro/^o = 0.427 provides 
the best overall fit to the observed wavelengths and wave speeds. On the adiabatically expanding 
jet the high initial sound speeds and relatively low initial Lorentz factor put L02 ~ uj* initially and 
ki{oJ2) > ki(uj\) for z b > 0.44 mas. 
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Fig. 12. — The upper panels indicate the growth rate, ki, of harmonic frequency pairs beginning at z b = 0.44 mas. 
The lower panels indicate the observed wavelength, \ ob , and wave speed, with the vertical lines indicating where 
z b ~ 25 mas. High frequency, LU2, (dotted line) and low frequency, uo\, (dashed line). 



Wave growth and damping for case B along with the accompanying pressure and velocity 
fluctuations is shown in Figure 13. The higher frequency with initial amplitude A2/R0 ~ 0.0035 
is damped as the more slowly growing low frequency wave is allowed to grow freely from its initial 
amplitude, Ai/Rq ~ 0.0013. Amplitude growth here should be compared to that shown in Figure 7 
for high frequency growth near to the core and isothermal expansion, and, in particular, to Figure 
9 for the development of a harmonic frequency pair on an isothermally expanding jet. Here initial 
amplitudes can be very small as the initial growth rates are high and 002 reaches saturation quickly. 
Amplitude growth is consistent with high frequency growth to saturation subsequently supplanted 
at larger distance by a more slowly growing lower frequency. Reduced velocity fluctuation at larger 
core distance compared to the isothermal model for maximal pressure fluctuation is a result of 
high Lorentz factor and low jet sound speed, i.e., high relativistic Mach number jMj. Overall, 
this dynamical result is more consistent with the required behavior of rapid high frequency growth 
followed by subsequent damping as the lower frequency grows to saturation. 

A pseudo-synchrotron intensity image for the adiabatically expanding jet is shown in Figure 
14. The apparent wave speed and wavelength undergoes a much more rapid initial increase than 
the isothermal Case 2 (Figure 8) from wavelength X oh < 1 mas at < 2 mas to X oh ~ 2.5 — 3.8 mas 
when z Q b ~ 2 — 8 mas. No additional wavelength or wave speed change for the higher frequency wave 
occurs beyond about 5 mas from the core. The low frequency wave dominates with X ob (uii) ~ 8 mas 
when z b > 14 mas. The apparent wavelength and wave speed of the adiabatic model fits the 
observed wavelength and motions in the inner 30 mas of the 3C 120 jet somewhat better than the 
isothermal model. However, here the intensity falls much more rapidly than is observed on the 
3C 120 jet. While neither model produces significant helical components at z Q b < 10 mas, the 
adiabatic jet shows no evidence for even modest southern helical components at z Q b > 15 mas 
even though the viewing angle is identical to the isothermal intensity image in Figure 10. The 
lack of helical components results from projection effects associated with relativistic aberration 
and a reduction in differential Doppler boosting. Here the helical pattern is viewed at an angle 
9 = 0.21 rad > 1//3 W when z Q b > 8 mas (see Figure 11). This has the effect of producing northern 
"helical" knots at 14 and 22 mas. This implies that the jet is at a viewing angle less than 12°. 
Differential Doppler boosting is reduced because the relativistic Mach number rises and transverse 
velocity fluctuation declines (see Figure 13). 
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Fig. 13. — The topmost panel of each set shows the growth or decline of A/Rj as a function of distance along the 
jet axis for a viewing angle of 9 = 12°. The dotted lines & dashed lines show A^IRj & Ai/Rj, respectively. Three 
panels below show the normalized pressure P/Po(z), and velocity components ^f z v z /c, v x /c and v y /c. The vertical 
dotted lines mark the projected separation from the core in mas. ID cuts are at x/Rj = 0.22 (dash), 0.44 (long 
dash), 0.66 (short-long dash), 0.88 (solid). 
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Fig. 14. — Pseudo-synchrotron intensity for an adiabatically expanding jet with harmonic frequencies at viewing 
angle 12° corresponding to the fluctuations shown in Figure 13. Contours are in factors of 2. 
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4.3. An Optimized Isothermal Image 

Constraints on our modeling are relaxed somewhat if much less acceleration in pattern speed 
is required, i.e., from about /3° < 4 when z Q b ~ 2 — 8 mas to about ~ 6 when z Q b > 14 mas. 
In this case the jet can be closer to isothermal than our present results suggest, although some 
sound speed decline and accompanying jet acceleration are still required. Such a result proves 
more consistent with the observed 5 GHz intensity decline along the jet which is about a factor 
of 3000 (see Figure 3). In order to make better comparison with the observed radio intensities we 
have constructed a maximal pseudo-synchrotron intensity image based on the isothermal Case 2 
harmonic frequency pair. This theoretical image along with images smoothed by Gaussian beams 
appropriate to 22 and 5 Ghz is shown in Figure 15. In constructing the theoretical intensity image 
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Fig. 15. — Pseudo-synchrotron intensity images for an optimized pressure fluctuation saturated isothermal expand- 
ing jet with harmonic frequencies at viewing angle 12°. The topmost is the initial computed theoretical image. The 
middle image is smoothed by a 0.3 x 0.6 mas Gaussian beam appropriate to 22 GHz, has a peak of 8126 units/beam 
and the lowest contour is 0.01 units/beam. The bottom is image smoothed by a 1.3 x 2.6 mas Gaussian beam appro- 
priate to 5 GHz, has a peak of 3949 units/beam, and the lowest contour is 0.1 units/beam. In all images contours 
are in factors of 2. The convolved images have considerably larger dynamic range than the observations. 

we have arbitrarily fixed the amplitude of the lower frequency at A\/Rj = 0.0035, i.e., achieves 
saturation at z b > 20 mas and allowed A2/R3 < A2/R0 ~ 0.04 to decline while maintaining 
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saturation of the harmonic combination. When z \y < 7 mas velocity and pressure fluctuations 
are similar to the Case 2 viewing angle 9 = 12° image shown in Figure 8. This maximizes the 
fluctuations at both small and larger scales. 

The 5 GHz image in Figure 15 shows a steeper intensity decline than seen in the observations. 
Some of this difference can be explained by core absorption at 5 GHz not included in the Figure 
15 image. Thus, the typical intensity decline as indicated by images at 5 GHz is on the order 
of isothermal. This result when combined with the somewhat better fit to overall dynamics and 
pattern speed provided by the adiabatic model suggests decoupling between macroscopic heating 
of the jet fluid as a whole, less than isothermal, and microscopic energization of the radiating 
electrons, at least isothermal. However, it is clear that adiabatic compressions are insufficient to 
give the observed helical components and the knot interknot variation observed (particularly at the 
higher 22 and 43 GHz frequencies). Different viewing angle can result in somewhat higher knot 
interknot intensity ratios at least in the innermost short wavelength region but not enough. Note 
that only viewing angles less than 12° are allowed as larger angles are ruled out by relativistic 
aberration effects. Thus, if helical components are entirely the result of projection and Doppler 
boosting they are not adiabatic compressions. The implication is that the leading edge of the helical 
pattern is a shock compression. Note that in this case the intensity might scale approximately with 
the jet density even if macroscopic heating is insufficient to make the jet fluid isothermal. We 
might also expect the intensity to be more ribbon like than the existing adiabatic compression as 
the pressure profile would be strongly peaked near to the jet surface. Furthermore, differential 
Doppler boosting effects may be enhanced by a flow field more aligned with the shock face, i.e., 
helical flow pitch more closely coupled to the intrinsic pitch of the helical shocked pressure and 
density enhancement. 

5. Summary & Discussion 

We have assumed that prominent components at z & ~ 2 — 30 mas are associated with a 
projected helical twist. These "helical" components indicate the wavelength and motion of the 
helical twist. Nearly constant spacing, X ob ~ 3 mas, and similar motion, ~ 4, of prominent 
components along with a "saturated" component intensity structure between 2 — 10 mas suggests 
that the motion of these components could be associated with a saturated helical wave in the high 
frequency regime. Between 10— 15 mas a discontinuity in the regular pattern of components occurs 
and between 15 — 30 mas the component spacing indicates a potential helical pattern with a twist 
wavelength X ob ~ 8 mas, or about twice the wavelength inside 10 mas. The increase in component 
spacing is accompanied by an apparent increase in the motion of components to /3° 6 ~ 6. Non- 
helical components outline a curved path in intensity contours between helical components and 
their motion is assumed indicative of the flow speed (3°^ ~ 6 — 6.5. 

In general wave growth, wavelengths, and wave motions are consistent with jet expansion 
that lies somewhere between adiabatic and isothermal but closer to adiabatic. Saturation can 
occur naturally on the expanding jet as a consequence of jet expansion and declining growth rate 
of helical waves at frequencies above resonance. Modeling the observed component spacing and 
motion as a helical twist from 2-30 mas requires two harmonic twist frequencies. The higher 
frequency has a period of < 2 years. A typical intrinsic wavelength of the superluminally moving 
helical twist produced by the two frequencies is six times the jet radius. Declining jet sound speed 
accompanied by modest jet acceleration provides the appropriate conditions for growth of multiple 
frequencies. Still it is necessary to assume a damping mechanism for the high frequency wave that 
allows the low frequency wave to grow to saturation. 

Modeling of the helical twist at the 6 = 12° viewing angle that has been extensively considered 
here, suggests a decline in jet sound speed from c/3 < ctj < cj\f?> at ~ 0.5 mas from the core to 
0.1c < dj < 0.25c at ~ 25 mas from the core accompanied by some acceleration in the jet flow from 
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Lorentz factor 7 < 5 to 7 > 7. At z \y ~ 15 mas /3 S = aj/c > 0.13 is obtained in the limit sin# — ► 0. 
The initial sound speed estimate at ~ 0.5 mas from the core increases to a,j ~ c/s/3 at a maximum 
possible viewing angle of 8 ~ 15°. At the larger viewing angle, higher typical Lorentz factor is 
required but additional jet acceleration also occurs. Thus, the initial Lorentz factor would not be 
too much higher than for the 12° viewing angle. However, viewing angles larger than 12° can be 
ruled out. At larger viewing angles relativistic aberration produces northern instead of southern 
"helical" knots when /3° 6 > 5 for > 10 mas. At smaller viewing angle the initial sound speed is 
less and less flow acceleration would occur. At a viewing angle of 6° the jet sound speed might be 
about half of the estimates at 12°, and the typical Lorentz factor would be comparable but with 
less variation along the jet. For z j, > 2 mas the flow Lorentz factor is in the range 7 ~ 6 — 8 for 
viewing angles 9 ~ 4 — 14°. 

Subluminal component motion at z {, < 1 mas provides an estimate of the sound speed in the 
cocoon medium outside the jet. The subluminal component motion requires an external (cocoon) 
sound speed a x > 0.1c at z Q b ~ 0.5 mas that is relatively independent of the viewing angle. The 
value depends on the Lorentz factor at this location. A lower Lorentz factor requires a higher 
external sound speed but less than the jet sound speed by a factor 1.5 - 5 depending on viewing 
angle. We cannot be more accurate as only the innermost z b < 1 mas of the jet provides any 
significant information on slow moving pattern speeds. This estimate was made by assuming that 
the innermost subluminal components were associated with the helical twist. However, similar 
results will be obtained if the components are modeled as trailing pinch structures triggered by 
passage of a shock down the jet (Agudo et al. 2001; Gomez et al. 2001). The combination of jet 
expansion, jet acceleration and sound speed decline implied by our results automatically guarantees 
that components associated with normal modes will show rapid acceleration and increased spacing 
at z f) < 2 mas. More modest increase in speeds and spacing occurs at larger core distance. 

We cannot determine how the external cocoon sound speed might vary with core distance 
from the observations considered here. Basically any determination of the change in external sound 
speed requires information on low frequency or resonant wave motion far from the core and here we 
see only high frequency waves (far above resonance) far from the core. Note that the present values 
for the external sound speed and an assumed adiabatic decline are consistent with a hot slow wind 
around a high speed hotter jet. Such a picture may be reasonable at the parsec to tens of parsec 
scales that we consider here. A wind with outflow speeds of a few tenths of lightspeed influences 
normal mode dynamics and increases the transverse velocity fluctuation produced by a saturated 
pressure fluctuation, e.g., Hardee & Hughes (2003). This effect would enhance differential Doppler 
boosting, possibly to significant levels. 

In earlier work Walker et al. (2001) showed that a flow helical pitch comparable to the intrinsic 
helical pitch implied by the observations can result in too much knot interknot intensity variation as 
a result of differential Doppler boosting. In this case differential Doppler boosting depended strongly 
on the viewing angle and thus knot interknot intensity variation could be used to constrain the 
viewing angle. Unfortunately, the present modeling of intensities as adiabatic compressions does 
not deliver sufficient helical projection plus differential Doppler boosting to produce the observed 
knot interknot intensity variations at any allowed viewing angle. Here our first order accurate flow 
field helical pitch is strongly decoupled and much less than the intrinsic helical pitch of the pressure 
and density enhancement. Differential Doppler boosting could be enhanced by a significant wind 
around the jet. Alternatively the apparent decoupling between the macroscopic fluid properties 
implied by an accelerating cooling jet and the microscopic radiating particle properties implied by 
an isothermal decline in the intensity can be taken to indicate the presence of a shock along the 
leading edge of the helix. The presence of a shock would more nearly align the flow field parallel to 
the shock face and enhance differential Doppler effects. The investigation of these possible effects 
is beyond the scope of this paper. 
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Finally our dynamical modeling implies that macroscopic viscosity and heating are relatively 
low along this jet, at least when compared to the energization requirements of the radiating electrons 
implied by our intensity modeling. A low macroscopic viscosity is consistent with the requirement 
that a relatively sharp velocity discontinuity between jet and cocoon must be present in order 
for high frequency short intrinsic wavelength helical patterns to exist. The requirement of sharp 
velocity discontinuity indicates a relatively low rate of mass entrainment. Broad velocity shear, 
indicative of higher macroscopic viscosity and mass entrainment, would eliminate coherent struc- 
tures based on short normal mode wavelengths, e.g., Birkinshaw (1991). Note however, that a 
shear layer increasing in thickness proportionally to the jet radius may be one reason for a growing 
lower frequency wave supplanting a higher frequency wave. The higher frequency shorter intrinsic 
wavelength would be damped as the jet expands and the shear layer grows in thickness while the 
lower frequency longer intrinsic wavelength could continue to grow. 

The European VLBI Network is a joint facility of European, Chinese, South African and other 
radio astronomy institutes funded by their national research councils. P. Hardee would like to thank 
Jean Eilek for numerous discussions concerning jet dynamics and acknowledges partial support 
from the NASA/NSSTC to the University of Alabama. J.L. Gomez acknowledges support from 
the Spanish Ministerio de Educacion y Ciencia and the European Fund for Regional Development 
through grant AYA2004-08067-C03-03. 

A. Observations at 5 GHz 

The 5 GHz images presented in Figure 2 are based on global VLBI observations of 3C 120 made 
on 1999 February 11 (1999.11) and on 1999 September 21 (1999.72). The data are from the ground 
stations involved in Halca orbiting VLBI observations. The data were recorded in 2 baseband 
channels of 16 MHz each using an array, in February, of 19 antennas; the Very Long Baseline Array 
(VLBA — 10 antennas), the phased Very Large Array, the Green Bank 140 foot, and 7 antennas of 
the European VLBI Network (Effelsberg, Jodrell Bank, Medicina, Noto, Onsala, Westerbork, and 
Sheshan). In September, the array was somewhat smaller with Green Bank, Onsala, Noto, and 
Sheshan not included. 

Correlation was done on the VLBA correlator and data reduction was done in AIPS. The gain 
normalization during self calibration was restricted to the VLBA antennas. The flux density scale 
is good to around 5%. A bandpass calibration was done to help achieve high dynamic range by 
minimizing closure errors. Many iterations of self calibration and imaging were required because 
of the complexity of the source — much of the structure is outside the windows displayed in this 
paper. The final off-source rms noise levels were 49 and 41 ^Jy beam^ 1 for the February and 
September images respectively. On-source errors, as usual, are expected to be significantly higher. 

We thank John Benson of NRAO for his significant contributions to the processing of these 
observations. The eventual main publication of these observations will show structures on signif- 
icantly larger scales than those displayed here. They also will include the much higher resolution 
images based on the Halca spacecraft data. 
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